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We report a direct measurement of D°-D° mixing parameters through a time-dependent amphtude 
analysis of the Dalitz plots of 7?" — ^ K°it'^-it~ and, for the first time, 7?" — ^ KgK'^K~ decays. The 
low-momentum pion tt^ in the decay Z)*+ — > D'^nf identifies the flavor of the neutral D meson 
at its production. Using 468.5 fb^^ of e'^e^ coUiding-beam data recorded near ^/s = 10.6 GeV 
by the BABAR detector at the PEP-II asyrnmctric-encrgy collider at SLAG, we measure the mixing 
parameters x = [1.6 ± 2.3 (stat.) ± 1.2 (syst.) ± 0.8 (model)] x 10"^, and y = [5.7 ± 2.0 (stat.) ± 
1.3 (syst.) ± 0.7 (model)] x 10~^. These results provide the best measurement to date of x and y. 
The knowledge of the value of x, in particular, is crucial for understanding the origin of mixing. 



PACS numbers: 13.25.Ft, ll.30.Er, 12.15.Ff, 14.40.Lb 

Particle-antiparticle mixing and CP violation (CPV) 
in the charm sector are predicted to be very small in the 
standard model (SM) [1 5]. Evidence for D"-D° mixing 
has been found only recently [6-11] and CPV has not 
been observed. Although precise SM predictions for D^- 
mixing are difficult to quantify, recent calculations of 



the mixing parameters x and y allow for values as large 
as O(10~^) [1]. The analyses to date that have reported 
evidence for mixing have not been able to provide direct 
measurements of x and y. A time-dependent amplitude 
analysis of the Dalitz plot (DP) of mesons decaying 
into K^-K^-K~ and K'^K^K~ self-conjugate final states 
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offers a unique way to access tlie mixing parameters x and 
y directly. In this Letter we study tlie time evolution of 
these three-body decays as a function of the position in 
the DP of squared invariant masses =771^(^^3/1+), s_ = 
rn?{K^h~), where h represents tt or K, and report the 
most precise single measurements of x and y to date. The 
knowledge of the value of .x, in particular, is crucial for 
understanding the origin of mixing and for determining 
whether contributions beyond the SM are present. 

We use the complete data sample of 468.5 fb"'^ 
recorded near ^/s = 10.6 GeV by the BaBAR experi- 
ment [12] at the PEP-II asymmetric-energy e+e" col- 
lider. The flavor of the neutral D meson at production 
is identified through the charge of the 7r+ ( "slow pion" ) 
produced in the decay D*+ 0°!^+ [13]. The £)° and 
mesons evolve and decay as a mixture of the Hamilto- 
nian eigenstates Di and D2, with masses and widths mi, 
Fi and m2, respectively. These mass eigenstates can 
then be written as linear combinations of flavor eigen- 
states, |L»i,2)=p|D°)±g|^°), where \p\^ + \q\'^ = 1. The 
mixing parameters are defined as a; = (mi — m2)/r and 
2/ = (Ti - r2)/2r, where T = (Pi r2)/2 is the average 
decay width. 

Assuming no CPV in the decay, the relation 
^(s_|_,s_) = ^(s_,s_|_) holds, where A and A arc the 
decay amplitudes for a D° or a D° into the final state 
K^h'^h~ as a function of the position in the DP. The 
time-dependent decay amplitude for a charm meson 
tagged at i = as D'^ or can then be written as 

M{s+, s-,t) = A{s+,s^)g+it) + ^Ais-,s+)g-{t), 
M{s+,s-,t) = fA{s+,s-)g+it) +A{s-,s+)g-{t), 

where g±{t) = 1/2 [e"*(™i-*ri/2)t _^ p-»(m2-»r2/2)tj 
q/p = 1 if CP is conserved in the mixing amplitude. 
The decay rates for and are obtained by squaring 
and A4 respectively, and consist of a sum of terms 
depending on (s+,s_) and proportional to cosh.{yTt), 
smh{yrt), cos(a;rt), and sin(a;rt), all modulated by the 
exponential decay factor e~^*. Assuming a model for 
A{s+,S-), it is possible to extract the mixing parame- 
ters X and y from the data, along with the amplitude 
model parameters and the proper-time resolution func- 
tion. The variation of the distribution of the events in the 
DP as a function of the proper D° decay time is the signa- 
ture of D'^-D^ mixing. The sensitivity to x and y arises 
mostly from regions in the DP where Cabibbo-favored 
and doubly-Cabibbo-suppressed amplitudes interfere and 
from regions populated by CP eigenstates [14]. This 
method was pioneered by CLEO [15] and extended to 
a significantly larger data sample by Belle [16]. 

The D° candidates are reconstructed in the if^Tr+Tr 
{KgK~^K~) final state by combining K'^ candidates with 
two oppositely-charged pions (kaons), with an invariant 
mass m^io between 1.824 and 1.904 GeV/c^. In order to 
reduce combinatorial background and to remove can- 



didates from i?-meson decays, we require the momen- 
tum of the D'-' in the e+e" center-of-mass frame to be 
greater than 2.5 GeV/c. The difference Am between the 
D*"*" and D'^ reconstructed invariant masses is required to 
satisfy 0.143 < Am < 0.149 GeV/c^. Each pion (kaon) 
track is identified using a likelihood particle identifica- 
tion algorithm based on cLE / dx ionization energy loss and 
Cherenkov angle measurements. The candidates are 
selected by pairing two oppositely-charged pions whose 
invariant mass is within 9 MeV/c^ of the nominal 
mass [17]. We require the cosine of the angle between 
the Kg flight direction (defined by the production 
and decay vertices) and the momentum to be greater 
than 0.99, and a decay length projected along the K'^ 
momentum to be greater than 10 times its error. These 
selection criteria suppress to a negligible level the back- 
ground from — >• ■K~^Tr~h~^h~ decays. For each charged 
track we require a transverse momentum with respect to 
the beam axis to be greater than 100 MeV/c, and for 
tracks from the decay we additionally require at least 
two hits in the two innermost layers of the silicon vertex 
tracker [12]. 

The D° proper time t, and its error at, are obtained 
through a kinematic fit of the entire decay chain which 
constrains the K'^ and pion (kaon) tracks to originate 
from a common vertex and also requires the D'^ and the 
7r+ candidates to originate from a common vertex, con- 
strained by the position and size of the e+e" interaction 
region. This reduces the contribution from — )■ K^K^ 
decays (affecting only KgTT~^'jT~) to 3% of the total back- 
ground. We retain candidates for which the proba- 
bility of the fit is greater than 0.01%, \t\ < 6 ps, and 
at < 1 ps. The most probable value for at is about 
0.2 (0.3) ps for /sTOtt+tt- {KIK+R-) signal candidates. 
For events where multiple D*^ candidates share one or 
more tracks, we keep the candidate with the highest 
probability. After applying all selection criteria, wc 
find 744 000 (96 000) KItt'+t:- {K°K+K-) candidates. 
Their m^jo and Am distributions are shown in Fig. 1 
and in [18]. 

The mixing parameters x and y are determined from 
an unbinned, extended maximum-likelihood fit to the 
if^TT+Tr" and K'^K~^K~ samples over the observables 
mjjo, Am, s+, S-, t, and at- First, the signal and back- 
ground yields are determined from a fit to m^o and Am 
distributions. For the subsequent fits, we restrict events 
to the signal region illustrated in Fig. 1, defined to lie 
within twice the measured resolution around the mean 
rriijo and Am values, and holding the signal and back- 
ground yields fixed to their signal window rescaled values. 
In the mixing fit, the rrijjo and Am shapes are excluded 
to minimize correlations with the rest of the observables. 
Our reference fit allows for mixing but assumes no CPV. 
We then allow for CPV as a cross-check of the mixing 
results. To avoid potential bias, the mixing results were 
examined only after the fitting and analysis procedures 
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FIG. 1: Distributions of m^o-Am for (a) KsTv^tt^ and (b) 
KsK~^K~ data after all selection criteria. The gray scale 
indicates the number of events per bin. The rectangles enclose 
the signal region. 



were finalized. 

For each fit stage, different sub-samples are character- 
ized separately: K^h'^h~ signal, random 7r+, misrecon- 
structed D°, ^ K°K° events (for ii'°7r+7r- decays 
only), and combinatorial background. The random 7r+ 
component describes correctly reconstructed D° decays 
combined with a random slow pion. Misreconstructed 
events have one or more decay products either missing 
or reconstructed with the wrong particle hypothesis. We 
account for K^K^ background events since they exhibit 
a characteristic DP distribution and a signal-like shape 
in the variables t and at- Combinatorial background 
events are those not described by the above components. 
The functional forms of the m/jo , Am probability density 
functions (PDFs) for the signal and background compo- 
nents are chosen based on studies performed on large 
Monte Carlo (MC) samples. These account for the ob- 
served correlations between Am and moo for signal and 
misreconstructed events. The PDF parameters are 
determined from two-dimensional likelihood fits to data 
over the full mjjo and Am region, or over dedicated side- 
band data samples. We find 540 800 ± 800 (79 900 ± 300) 
signal events in the if^Tr+Tr^ {K^K^ K~) signal region, 
with purities of 98.5% (99.2%), and reconstruction effi- 
ciencies of 14.4% (14.6%). Random 7r+, misreconstructed 
, and combinatorial background events account for 
23% (53%), 52% (23%), and 22% (24%) of the total back- 
ground. Projections of these fits and the contributions 
from the different background components can be found 
in [18]. 

The amplitudes A{s+, s_) are described by a coherent 
sum of quasi- two-body amplitudes [14, 19]. The dynam- 
ical properties of the P- and D-wave amplitudes are pa- 
rameterized through intermediate resonances with mass- 
dependent rclativistic Breit-Wigner (BW) or Gounaris- 
Sakurai (GS) propagators, Blatt-Weisskopf centrifugal 
barrier factors, and Zemach tensors for the angular dis- 
tribution [19]. The TTTT S-wave dynamics is described 
through a K-matrix formalism with the P-vector ap- 



proximation and 5 poles [14, 20]. For the Kir S-wave 
we include a BW for the /sTg (1430)=*= with a coherent 
non-resonant contribution parametrized by a scattering 
length and effective range similar to those used to de- 
scribe K-K scattering data [18, 21]. For the KK S-wave, 
a coupled-channel BW is used for the ao(980) isovector 
with BWs for the /o(1370) and ao(1450) states. 

We define PDFs to describe the dependence of the com- 
ponents in our event sample upon DP position (s_|_,s_) 

and upon decay time, t. For signal, \M\^ ov \m\ is 
convolved with a proper-time resolution function differ- 
ent for /ST^TT+Tr" and K^K^K~ events with parameters 
determined by the mixing fit to the data. The resolu- 
tion function is a sum of three Gaussians with one of the 
means allowed to differ from zero (the offset, to), and 
two of the widths proportional to at- While to does not 
depend sensitively on the DP position, the ability to re- 
construct t varies as a function of (s+,s_). Hence, the 
observed distributions of at in a number of DP regions 
arc included in the signal PDF. We apply corrections for 
efficiency variations and neglect the invariant mass res- 
olution across the DP. The time-dependent PDF for the 
small random tt^ background component is described by 
an equal combination of and £)° signal events assum- 
ing no mixing, since the slow pion is positive or neg- 
ative with approximately equal probability, and carries 
little weight in the vertex fit. The PDFs for misrecon- 
structed events and combinatorial background are de- 
termined from m]ja sideband samples. A non-parametric 
approach is used to construct the DP distributions, while 
the proper-time distributions are described by a sum of 
two Gaussian functions, one of which has a power-law tail 
to account for a small long-lived component. The back- 
ground components containing real and misreconstructed 
decays have different at distributions, which are de- 
termined from the signal and mj^o sideband regions. 

Results for our nominal mixing fit, in which D'^ and 
samples from /C^Tr+Tr" and K^K^K~ channels are 
combined, are reported in Table I. The proper-time dis- 
tributions with their fit projections are shown in Fig. 2. 
Additional fit results and projections can be found in [18]. 
We evaluate the amplitude model fit to the DP distribu- 
tion with a test with two-dimensional adaptive bin- 
ning, and obtain = 10 429.2 (1 511.2) for 8 626 - 41 
(1 195 — 17) degrees of freedom (ndof), for K^tt^'it~ 
[K^K^ K~). MC studies show that a significant contri- 
bution to these x^ values, Ax^/ndof ~ 0.16, arises from 
imperfections in modeling experimental effects, mostly 
the efficiency variations at the boundaries of the DP and 
the invariant mass resolution. The fitted average lifetime 
T = l/F is found to be consistent with the world av- 
erage lifetime [17], while to is found to be 5.1 ± 0.8 fs 
(5.1 ±2.2 fs) for the /r°7r+7r- [K^K+K-) mode, consis- 
tent with expectations from small misalignments in the 
detector [10]. 
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FIG. 2: (color online) Proper-time distributions for (a) 
Ks-K^-K~ and (b) K^K'^ K~ data in the signal region with 
— 2 < t < 4 ps (points). The curves show the fit projec- 
tions for signal plus background (solid lines) and for different 
background components (shaded regions) . 



TABLE I: Results from the mixing fits. The first uncertainty 
is statistical, the second systematic and the third systematic 
from the amplitude model. For the nominal fit, the corre- 
sponding correlation coefficients between x and y are 3.5%, 
f6.0% and -2.7%, respectively. 

Fit type x/lQ'^ y/10~^ 



Nominal 
KlK+K- 



1.6±2.3±i.2±0.8 
2.6 ±2.4 
-13.6 ±9.2 
0.0 ±3.3 
3.3 ±3.3 



5.7 ±2.0 ±1.3 ±0.7 
6.0 ±2.1 

4.4 ± 5.7 

5.5 ±2.8 
5.9 ±2.8 



A variety of studies using large MC samples with both 
parameterized and full detector simulations and data 
have been performed to validate the analysis method 
and fitting procedure and to check the consistency of 
the results. These studies demonstrate that the analy- 
sis correctly determines the mixing parameters with in- 
significant biases and well-behaved Gaussian errors. No 
significant variations of the mixing parameters are ob- 
served as a function of momentum, polar and azimuthal 
angles of the 13° meson, and data taking period. In- 
cluding the m£)0, Am PDFs in the mixing fit does not 
significantly change the values for x and y. The mixing 
fit has also been performed separately for the X^Tr+Tr^ 
and K^K'^K^ data samples, and for and decays, 
with the results listed in Table I. Fitting separately for 
and provides a check against possible effects from 
CPV in mixing and in decay. Finally, if we fit the data 
forcing the decay amplitudes for and to be the 
same (no direct CPV), but allowing their x and y values 
to differ, we find these values to be consistent (i.e., no 
evidence for CPV in mixing). 

Systematic uncertainties arise from approximations in 
the modeling of experimental and selection criteria ef- 
fects [18]. We account for variations in the signal and 
background yields, the efficiency variations across the 



DP, the modeling of the DP and proper-time distribu- 
tions for events containing misreconstructed decays, 
the misidentification of the flavor for signal and ran- 
dom 7r+ events, potential effects due to mixing in the 
random 7r+ background component, and PDF normal- 
ization. We also consider variations of the resolution 
function and at PDFs, including alternatives to describe 
the correlation between at and the DP position (e.g. ne- 
glecting the dependence of the at distributions on the DP 
position entirely). The dominant sources of experimen- 
tal systematic uncertainty are the limited statistics of full 
detector simulations (used to study potential biases due 
to the event selection, invariant mass resolution, residual 
correlations between the fit variables, and fitting pro- 
cedure) and instrumental effects arising from the small 
misalignment of the detector. Effects from our selection 
criteria are estimated by varying the m^jo, Am, t, and 
at requirements. 

Assumptions in the amplitude models are also a source 
of systematic uncertainty [14, 18]. We use alternative 
models where the BW parameters are varied accord- 
ing to their uncertainties or changed to values mea- 
sured by other experiments, the reference K-matrix so- 
lution [14] is replaced by other solutions [20], and the 
standard parameterizations are substituted by other re- 
lated choices. These include replacing the GS by BW 
lineshapes, removing the mass dependence in the P- 
vector [22], changes in form factors such as variations 
in the Blatt-Weisskopf radius and the effect of evaluat- 
ing the momentum of the spectator particle in the 
meson frame rather than in the resonance rest frame, 
and adopting a helicity formalism [19] to describe the 
angular dependence. Other models are built by remov- 
ing or adding resonances with small or negligible frac- 
tions. The largest effect arises when the uncertainties in 
the amplitude model parameters obtained from the fit to 
the DP variables only [18] are propagated to the mixing 
fit. These uncertainties are dominated by the parameters 
related to the Kit S and P waves. 

The mixing significance is evaluated by the variation of 
the negative log-likelihood (— 2Aln£) in the mixing pa- 
rameter space. We account for the systematic uncertain- 
ties by approximating £ as a two-dimensional Gaussian 
with covariance matrix resulting from the sum of the cor- 
responding statistical, systematic, and amplitude model 
matrices. Figure 3 shows the confidence-level (C.L.) con- 
tours in two dimensions (x and y) with systematic un- 
certainties included. The variation in — 2Aln£ for the 
no-mixing point is 5.6 units which corresponds to a C.L. 
equivalent to 1.9 standard deviations, including the sys- 
tematic uncertainties. 

In summary, we have directly measured the mix- 
ing parameters x — [1.6 ± 2.3 (stat.) ± 1.2 (syst.) ± 
0.8 (model)] x lO^^, and y ^ [5.7 ± 2.0 (stat.) ± 
1.3 (syst.) ±0.7 (model)] x 10"'^, using, for the first time, 
a combined analysis of — > if^Tr+Tr" and — 



7 




-10 10 



FIG. 3: Central value (point) and C.L. contours (including 

statistical, systematic and amplitude model uncertainties) in 
the x-y plane for C.L. = 68.3%, 95.4%, 99.7%. The no-mixing 
point is shown as a plus sign (+). 



K^K^K~^ decays. These results are consistent with the 
previous similar measurements of -fCgTr+Tr" alone [15, 16] 
and have improved precision. They disfavor the no- 
mixing hypothesis with a C.L. equivalent to 1.9 stan- 
dard deviations and are in agreement with the range 
of SM predictions [1- 5]. Our measurements favor lower 
values for x than for and lower x and y values (but 
still consistent) than those obtained when combining re- 
sults from other decays [6 -11, 23]. Adding our re- 
sults to the combination of all previous analyses signif- 
icantly improves our current knowledge of the mixing 
parameters x and y, whose average values change from 
(9.8±2.5) X 10-3 and (8.3±1.6) x IQ-^ to (5.9±2.0) x lO'^ 
and (8.0 ± 1.3) x 10"^, respectively [24]. 
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TABLE I: D° 



""tt complex amplitudes, -ktt P-vector and Kir S-wave parameters, and fit fractions, as obtained from the 



mixing fit. The irn S-wave parameters ps, ffl""^, and ffr°'^ are fixed to zero due to the lack of sensitivity. We also report the 
mass and the width of the K* (892)^ resonance. Errors arc statistical only. The fit fraction is defined as the integral over the 
entire DP of a single component divided by the coherent sum of all components. The sum of fit fractions is 103.3%. A detailed 
description of the parameters can be found elsewhere [14]. Equations (14) and (15) in [14] have been corrected as follows 

Aktt L=o{s) = Tktt l 



o(s)/p(s), where p(s) = g/\/s is the phase-space factor and Tktt i,=o(s) = Fsm{SF + (f>F)e^^^^'^'^'''' + 



ilsin5fle'(*«+*«)e'2(«^+*^), with tanfe = MK.(i430)rK-(i430)(s)/(M| 



o(1430) 



- s), cot 5f = l/{aq)+rq/2, s the invariant mass 



squared of the Kn system, and q the momentum of the kaon (or pion) in the Kiv rest frame [21]. The symbol ^ indicates the 
parameters fixed in the mixing fit to the values extracted from a time-integrated DP fit to the same data. The results from 
this time-integrated DP fit for the amplitude model parameters agree within statistical errors with the results reported here. 



Component 


Amplitude 


Phase (rad) 


Fit fraction (%) 


K*{892)- 


1.735 ± 0.005 


2.331 


± 


0. 


.004 


57, 


.0 


p(770)« 


1 











21, 


.1 




2.650 ±0.015 


1.497 


± 


0, 


.007 


6, 


.1 




1.303 ± 0.013 


2.498 


± 


0. 


.012 


1, 


,9 


w(782) 


0.0420 ± 0.0006 


2.046 


± 


0. 


.014 


0, 


,6 


K*{892)+ 


0.164 ±0.003 


-0.768 


± 


0, 


,019 


0, 


,6 


A:*(1680)- 


0.90 ±0.03 


-2.97 


± 


0, 


,04 


0, 


,3 


/2(1270) 


0.410 ±0.013 


2.88 


± 


0, 


,03 


0, 


,3 


K5{U30)+ 


0.145 ± 0.014 


1.78 


± 


0. 


,10 


< 


.1 




0.115 ±0.013 


2.69 


± 


0. 


,11 


< 


.1 


\ nn S-wave 












15, 


,4 


/3i 


5.54 ±0.06 


-0.054 


± 


0, 


,007 






ft 


15.64 ± 0.06 


-3.125 


± 


0, 


,005 






ft 


44.6 ± 1.2 


2.731 


± 


0, 


,015 






ft 


9.3 ±0.2 


2.30 


± 


0. 


,02 






I'prod 

111 


11.43 ±0.11+ 


-0.005 


± 


0, 


.009+ 






/■prod 
Jl2 


15.5±0.4t 


-1.13 


± 


0, 


.02+ 






/■prod 
/l3 


7.0 ± 0.7+ 


0.99 


± 


0, 


.11+ 







prod 
®0 



Parameter value 
-3.92637 



[ Kn S-wave parameters 
Mk* (14.30) (MeV/c") 
rx*(i43(,) (MeV/c^) 
F 

4>F (rad) 
R 

(f)R (rad) 
a (GeV/c"') 
[r (GeV/c"') 



1421.5 ± 1.6+ 
247 ± 3+ 
0.62 ±0.04+ 
-0.100 ±0.010+ 
1 

1.10 ±0.02+ 
0.224 ± 0.003+ 
- 15.01 ±0.13+ 



[ K* (892) parameters 

Mk* (892) (MeV/c^) 
LrK*(892) (MeV/c^) 



893.70 ± 0.07+ 
46.74 ±0.15+ 
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TABLE II: D K^K^ K~ complex amplitudes and fit fractions, as obtained from the mixing fit. We also report the mass 
and the width of the cj>{lQ2Q) resonance, and the ao(980) coupling constant to KK as determined from the fit. Errors are 
statistical only. The fit fraction is defined as the integral over the entire DP of a single component divided by the coherent sum 
of all components. The sum of fit fractions is 163.4%. A detailed description of the parameters can be found elsewhere [14]. 
The symbol ^ indicates the parameters fixed in the mixing fit to the values extracted from a time-integrated DP fit to the same 
data. The results from this time-integrated DP fit for the amplitude model parameters agree within statistical errors with the 
results reported here. 



Component 


Amplitude 


Phase (rad) 


Fit fraction (%) 


ao(980)" 


1 





51.8 


0(1020) 


0.2313 ±0.0011 


-0.977 ±0.008 


44.1 


a,, (1450) + 


0.93 ±0.03^ 


1.66 ± 0.07+ 


25.6 


a„(980)+ 


0.635 ± 0.006 


-2.91 ±0.02 


19.5 


a„(1450)° 


0.83 ±0.10+ 


-1.93 ±0.12+ 


19.3 


/„(1370) 


0.16 ±0.05+ 


0.2 ±0.2+ 


1.7 


/2(1270) 


0.385 ±0.015 


0.06 ±0.04 


0.7 


ao(980)~ 


0.125 ±0.008 


2.47 ±0.04 


0.7 



[ (/>(1020) and ao(980) parameters Value 

Af0(io2O) (MeV/c^) 1019.55 ±0.02+ 

r,^(i020) (MeV/c^) 4.60 ±0.04+ 

Lfeif(MeV/c2) 537 ±9+ 




FIG. 1: (color online) Distributions of m^o and Am for (a,b) K'^-k'^-k~ and (c,d) K^K^ K~ data after all selection criteria 
(points). The curves superimposed represent the fit projections for signal plus background (solid lines) and for different 
background components (shaded regions). The arrows indicate the definition of the signal region. 
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s. (Gey^lc^) (GeV^/c*) 

FIG. 2: DP distributions for (a) D° -> Tv'^n and (b) D" — ^ K^K^ K data after all selection criteria, in the signal region. 
The gray scale indicates the number of events per bin. The solid lines show the kinematic limits of the decay. The so DP 
variable is defined as so = iv?{K^h~). For decays the variables s_ and s+ are interchanged. 




s. (GeV^/c*) s. (GeV^/c") So (GeV^/c") 



FIG. 3: DP projections for (a,b,c) D — >■ KsTv^tt and (d,e,f) D — > KsK^ K data after all selection criteria, in the signal 
region (points). The histograms represent the mixing fit projections. For decays the variables s_ and s+ are interchanged. 



TABLE III: Summary of the contributions to the experimental systematic uncertainty on the mixing parameters. 



Source 




y/10-3 


Analysis biases and fitting procedure (Monte Carlo statistics) 


0.75 


0.66 


Selection criteria 


0.47 


0.57 


Signal and background yields 


0.11 


0.07 


Efficiency variations across the DP 


0.37 


0.18 


Modeling of the DP distributions for misreconstructed decays 


0.33 


0.14 


Modeling of the proper-time distributions for signal and misreconstructed decays 


0.13 


0.13 


Modeling of the proper-time error distributions for signal and misreconstructed decays 


0.06 


0.09 


Misidentification of the flavor for signal and random tt^ events 


0.49 


0.40 


Mixing in the random irt background component 


0.10 


0.08 


PDF normalization 


0.11 


0.05 


Misalignment of the detector 


0.28 


0.83 


Total experimental systematic uncertainty 


1.18 


1.30 
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TABLE IV: Summary of the contributions to the D° decay amphtude model systematic uncertainty on the mixing parameters. 



Source 




J//10-3 


Breit-Wigncr parameters and alternative GS lineshapes 


U.oo 


U.iz 


Alternative K-matrix solutions and P-vector parameterization 


0.13 


0.19 


Ktv S- and P-waves, and tttt S-wave parameters 


0.68 


0.53 


Form factors 


0.25 


0.23 


Angular dependence 


0.05 


0.17 


Add/remove resonances 


0.17 


0.23 


Total amplitude model systematic uncertainty 


0.83 


0.69 



